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Results are reported based on samples of 58 million J/ip and 14 million ip(2S) de- 
cays obtained by the BESII experiment. Improved branching fraction measurements are 
determined, including branching fractions for J/ip — » 7r + 7r _ 7r°, ip(2S) — > ir°J/ip, r)J/ip, 
7r°7r° J/ip, anything J/ip, and ip(2S) — > 7Xcii7Xc2 — * "fyJ/^P- The decay J/ip — > u)7r+7r — 
is studied. At low -kit mass, a large, broad peak due to the a is observed, and its pole posi- 
tion is determined. Results are presented on ip(2S) and J/ip hadronic decays to KgpK~n 
and KgpK + n final states. No significant ©(1540) signal, the pentaquark candidate, is 
observed, and upper limits are set. An enhancement near the m p + M\ mass threshold is 
observed in the pA invariant mass spectrum from J/if) — » pK~ A+c.c. decays. It can be fit 
with an S-wave Breit-Wigner resonance with a mass m = 2075 ± 12 (stat) ± 5 (syst) MeV 
and a width of T = 90 ± 35 (stat) ± 9 (syst) MeV. 

Keywords: charmonium; pentaquark; hadronic transitions. 



1. Introduction 

The Beijing Spectrometer (BES) is a general purpose solenoidal detector at the 
Beijing Electron Positron Collider (BEPC). BEPC operates in the center of mass 
energy range from 2 to 5 GeV with a luminosity at the J/ip energy of approximately 
5 x 10 30 cm^s" 1 . BES (BESI) is described in detail in Ref.ffl and the upgraded 
BES detector (BESII) is described in Ref. H 



2. B(J/ip -> 7r+7r-7r ) 

The largest J ftp decay involving hadronic resonances is J/ip — > p(770)7r. Its branch- 
ing fraction has been reported by many experimental groups^ assuming all 7r + 7r~7r° 
final states come from p{710)ir. The precision of these measurements varies from 
13% to 25%. Here, we present two independent measurements of this branching 
fraction using J/ip and ip(2S) decays. 



1 



February 7, 2008 7:36 WSPC/INSTRUCTION FILE meson04'bes'fah 



2 Frederick A. Harris 

2.1. Absolute measurement of J/ip — ► 7r+7r _ 7r° decays 

Events with two oppositely charged tracks and at least two good photons are se- 
lected. A 5-constraint (5C) kinematic fit is made under the ir + ir~j-f hypothesis 
with the invariant mass of the two photons being constrained to the tt° mass, and 
the fit x^+jt-tt 1S rec l u i re d to be less than 15. After these requirements, 219691 
7T + 7r~7r° candidates are selected. The branching fraction is B(J/ip — > 7r + 7r~7r°) = 
(21.84 ± 0.05 ± 2.01) x 10~ 3 , where the first error is statistical and the second 
systematic. 

The Dalitz plot of 771^+^0 versus m^^o is shown in Fig. ^ Three bands are 
clearly visible in the plot, corresponding to J/ip — > pir; J/ip — > 7r + 7r~7r is strongly 
dominated by pn. 




5 4 6 

m 2 (7irc°) (GeV/c 2 ) 2 



Fig. 1. The Dalitz plot for J/ip -> tt+tt^tt' 



2.2. Relative measurement of J/ip — > 7T+7T 7r° 

The relative measurement is based on a sample of 14 million ip{2S) events. The 
ip(2S) is a copious source of J/ip decays: the branching fraction of ip(2S) — » 
J/ip is the largest single ip(2S) decay channel. Therefore, we can determine 
the branching fraction of J /ip — ► 7r + 7r~7r° from a comparison of the following two 
processes: 

ip(2S) -> 7T+7T- J/tP 

^ TT+TT-TT (I) 

and <—> /J + /i~ (II) 

Using the relative measurement, many systematic errors mostly cancel. Therefore, 
the precision of the branching fraction J/ip — > 7r + 7r~7r from the relative measure- 
ment is comparable with that of the direct J/ip decay, although the size of ip(2S) 
sample is smaller. 

For process I, a 5C kinematic fit is performed for each candidate event, and the 
event probability given by the fit must be greater than 0.01. A 4C kinematic fit is 
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performed for ip(2S) — ► Tt + tt~ (i + fi~ candidate events, and the probability given by 
the fit must be greater than 0.01. The branching fraction is B(J/ip — > 7r + 7r~7r°) = 
(20.91 ±0.21 ±1.16) x 10~ 3 . 

The results of the two measurements are in good agreement. Their weighted 
mean is 

B(J/t/j -► tt+tt-tt ) = (2.10 ± 0.12)%. 

The result obtained is higher than those of previous measurements and has better 
precision. For more detail, see Ref. 0] 

3. J/tp — ► u>ir+ir~ 

There has been evidence for a low mass pole in the early DM2 - and BESI 3 data on 
J/t/j — > um + ir~ . Here results on J/tp — > aj7r + 7r~ from 5.8 x 10 7 J/ip events collected 
with the BESII detector are presented. 

The u> is observed in its 7r + 7r~7r° decay mode. Events are required to have four 
good charged tracks with total charge zero and more than one good photon. The 
TOF and dE/dx information are used to identify pions; they largely reject kaons 
from background reactions such as K + K~ 7r + ir~7r°. Events with a 27 invariant mass 
M 77 - M w o\ < 40 MeV/c 2 are fitted with a 5C kinematic fit to tt + tt^tt + tt^tt with 
the two photons being constrained to the ir° mass. Events with xtc < ^0 are 
selected. The resulting 7r + 7r~7r° mass distribution for is shown in Fig. Et a )- The ui 
signal is selected by requiring [M^+^-^o — M u \ < 40 MeV/c 2 . 

Fig-Etb) shows the ir + ir~ invariant mass spectrum which recoils against the u>, 
and Fig. Etc) shows the uiir invariant mass. The Dalitz plot of this channel is shown 
in Fig.|2fd). There is a large ^2(1270) peak in Fig.Etb) and a strong &i(1235) peak 
in Fig. Etc). At low tttt masses in Fig. Etb), a broad enhancement which is due to 
the a pole is clearly seen. This peak is evident as a strong band along the upper 
right-hand edge of the Dalitz plot in Fig.Etd). 

Partial wave analyses (PWA) are performed on this channel using two methods. 
In the first method, the whole mass region of M w +^- which recoils against the u> 
is analyzed, the u> decay information is used, and the background is subtracted 
by sideband estimation. For the second method, the region M 7r + lr - < 1.5 GeV is 
analyzed, and the background is fitted by bn phase space. In both methods, different 
parameterizations of the a pole are also studied. 

The mass and width of the a are different when using different a parameteriza- 
tions. However, the pole position of the a is stable; different analysis methods and 
different parameterizations of the a amplitude give consistent results for the a pole. 
From a simple mean of the six analyses, the pole position of the a is determined to 
be (541±39 - i (252±42)) MeV. Here, the errors cover the statistical and systematic 
errors in the six analyses, as well as the error in the extrapolation to the pole. The 
systematic errors dominate. More detail may be found in Ref. [7\ 
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4. i/)(2S) -> -y-yJ/ip and tp{2S) ->• XJ/i/j 

Experimental results for the processes V^S) ~ > ^J/ip, wJ/ip, and 7Xci,2 are few 
and were mainly taken in the 1970s and 80s. ^ Here, we report on the analysis of 
i(j{2S) — » ir°J/ip, T/J/ip, and 7X c i,2 decays based on a sample of 14.0 x 10 6 0(25*) 
events collected with the BESII detector. Events with two charged tracks identified 
as an electron pair or muon pair and two or three photon candidates are selected. 
A five constraint (5C) kinematic fit to the hypothesis '0(25') — * 77? + Z - with the 
invariant mass of the lepton pair constrained to J/ip mass is performed, and the fit 
probability is required to be greater than 0.01. 

Fig. El shows, after a cut to remove the huge background from tp(2S) —> jXci,c2 
under the ip(2S) — > ir°J/ip signal, the distribution of invariant mass, M 77 . A Breit 
Wigner with a double Gaussian mass resolution function to describe the ir° reso- 
nance plus a background polynomial is fitted to the data. Similar analyses are made 
for the other channels. 

Using the fitting results and the efficiencies and correction factors for each chan- 
nel, the branching fractions listed in Tableware determined. The BES B(ip(2S) — > 
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Fig. 3. Two photon invariant mass distribution for candidate p(2S) — > n J/tp events for (a) 
77e~'~e _ and (b) -yyfi + ii~. 





Table 1. Results from ip(2S) - 


» nJ/4> 


Channel 


TT U J/lp 




Final state 


77e + e — 77/j+/j~ 


77e + e _ 77^ + M~ 


B (%) 
Combined (%) 
PDG (%P 


0.139 ± 0.020 ± 0.013 0.147 ± 0.019 ± 0.013 
0.143 ±0.014 ±0.013 
0.096 ± 0.021 


2.91 ±0.12 ±0.21 3.06 ±0.14 ±0.25 
2.98 ±0.09 ±0.23 
3.16 ±0.22 


Channel 


7Xci 


7Xc2 


Final state 


77e + e — 77 fi + fi~ 


77e + e _ 77/i + /^ — 


B(%) 
Combined (%) 
PDG (%P 


8.73 ±0.21 ± 1.00 9.11 ±0.24 ± 1.12 
8.90 ±0.16 ±1.05 
8.4 ± 0.6 


7.90 ±0.26 ±0.88 8.12 ± 0.23 ± 0.99 
8.02 ±0.17 ±0.94 
6.4 ±0.6 



TT°J/%li) measurement has improved precision by more than a factor of two com- 
pared with other experiments, and the ip(2S) — > r/J/ip branching fraction is the 
most accurate single measurement. More details on this analysis may be found in 
Ref.El 

In another analysis, based on a sample of approximately 4 x 10 6 ip{2S) events 
obtained with the BESI detector, ^ a different technique is used for measuring 
branching fractions for the inclusive decay ip{2S) — > anything,//?/), and the exclusive 
processes for the cases where X = 77 and X = irir, as well as the cascade processes 
ip(2S) — > 7XeO/i/2, XeO/1/2 — ► jJ/ip- Inclusive pairs are reconstructed, and 

the number of rp(2S) — > J/tJjX events is determined from the J/ip — > peak 
in the invariant mass distribution. The exclusive branching fractions are 

determined from fits to the distribution of masses recoiling from the J '/tp with 
Monte-Carlo determined distributions for each individual channel. 

The mass recoiling against the J/ip candidates, mx is determined from energy 
and momentum conservation. To determine the number of exclusive decays and 
separate ip(2S) — > J/i/^tt ^ and ip(2S) — > J/-07r + 7r~ events, mx histograms for 
events with and without additional charged tracks, shown in Figs. 0] and [5J are fit 
simultaneously. To reduce background and improve the quality of the track momen- 
tum measurements, events used for this part of the analysis are required to have 
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a kinematic fit % 2 < 7. The channels of interest are normalized to the observed 
number of 7r + 7r~ J/ip events; ratios of the studied branching fractions to that for 
B(ip(2S) — > 7r + 7r - J/ip) are reported. This has that advantage that many of the 
systematic errors largely cancel. 




Fig. 4. Fit of the rax distribution events 
with no additional charged tracks. Shown 
are the data (points with error bars), the 
component histograms, and the final fit. For 
the components, the large, long-dash his- 
togram is ^}{2S) — * J/ipirn, the narrow, 
dash-dot histogram is ip(2S) — > J/ipVi the 
broad, short-dashed histogram is ip(2S) — » 
7Xcl>Xcl - * ~fJ/ip, the broad, hatched his- 
togram is i))(2S) —> 7Xc2,Xc2 — * jJ/tp, and the 
lowest cross-hatched histogram is the combined 
e + e - _> and e+e~ -+ i/)(2S), ip(2S) -» 

(■y)fi + H~ background. The final fit is the solid 
histogram. 



Fig. 5. Fit of the rax distribution for events 
with any number of additional charged tracks. 
Shown are the data (points with error bars), 
the component histograms, and the final fit 
(solid histogram). The dashed histogram is 
ip{2S) — » J/»/)7r+7r _ , and the hatched his- 
togram is ip(2S) — > J/ipr/. There is very little 
evidence for ip(2S) -» 7Xci/2>X c i/2 -* "fJ/f- 
This distribution is composed predominantly of 
ip(2S) -> J/tpn+TT-. 



The final branching fraction ratios and branching fractions are shown in Tabled 
along with the PDG experimental averages and global fit results. The results for 
B(J/ip anything) /B(ip(2S) -» n+w-J/ip) and B(r]J/ip)/B(tp(2S) -» tt+tt-J/^) 
have smaller errors than the previous results. The agreement for both the ratios of 
branching fractions and the calculated branching fractions using the PDG result for 
ip(2S) — > B(tt + tt~ J ftp) with the PDG fit results is good, and the determination of 
B(r)J/ip) agrees well with the determination from rp(2S) — > ■yfJ/tp decays above. 
More details on this analysis may be found in Ref. 
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Table 2. Final branching ratios and branching fractions. PDG04-exp results are single 
measurements or averages of measurements, while PDG04-fit are results of their global fit 
to many experimental measurements. The BES results in the second half of the table are 
calculated using the PDG04 value of B nw = B(i>(2S) -» J/i/)7r+7r _ ) = (31.7 ± 1.1)%. 



Case 


This result 


PDG04-cxp 


PDG04-fit 


B( J/ip (anything)/B^^ 
B(J/V7r°7r°)/B w ^ 
BiJ/^/B^ 
B(7Xd)B(Xci - 1 J /ip)/ 
B(rXc2)B( Xc2 -* jJ/^)/B nn 


1.867 ±0.026 ±0.055 
0.570 ± 0.009 ± 0.026 
0.098 ±0.005 ±0.010 
0.126 ±0.003 ±0.038 
0.060 ± 0.000 ± 0.028 


2.016 ±0.150 

0.091 ±0.021 
0.085 ±0.021 
0.039 ±0.012 


1.821 ±0.036 
0.59 ±0.05 
0.100 ±0.008 
0.084 ± 0.006 
0.041 ± 0.003 


B{J/iP (anything) (%) 
B(J/tPtt°tt ) (%) 
B(J/rh) (%) 

B(7Xc l )B(Xci-»'YJ r /V0 (%) 
S(7Xc2)B(Xc2 (%) 


59.2 ±0.8 ±2.7 
18.1 ±0.3 ±1.0 

3.11 ±0.17 ±0.31 
4.0 ±0.1 ±1.2 

1.91 ±0.01 ±0.86 


55 ±7 

2.9 ±0.5 
2.66 ±0.16 
1.20 ±0.13 


57.6 ± 2.0 
18.8 ± 1.2 
3.16 ±0.22 
2.67 ±0.15 
1.30 ±0.08 



5. Pentaquark Search 

The pentaquar the main topic of the opening session of this conference, has 
generated much excitement. BES has searched for the pentaquark state 0(1540) in 
i>(2S) and J /ip decays to KgpK~n and KgpK + n final states using samples of 14 
million ip(2S) and 58 million J/ip events taken with BES II. These processes could 
contain O decays to Kgp, K + n (uudds) and decays to Kgp, K~n (uudds). 
The anti-neutron and neutron are not detected. The Kg meson in the event 



is identified through the decay Kg 



Candidate events are kinematically 



fitted under the assumption of a missing h(n) to obtain better mass resolution and 
to suppress the backgrounds. Events with missing mass close to the n(n)'s mass are 
selected. We use the same criteria and treatment for both ip(2S) and J/ip data. 



> 
o 



1.8 2 2.2 

M(K s p(K + n)) (GeV/c 2 ) 



Fig. 6. Scatter plot of K~n {K° s p) versus K°p (K+n) for ip(2S) 
modes. 



K%pK-n + K°pK+r 



The scatter plot of K~n (K° s p) versus K° s p (K+n) for ^(25) -» K%pK-fi + 
KgpK + n modes is shown in Fig. EI Zero events fall within the signal region, shown 



7, 2008 7:36 WSPC/INSTRUCTION FILE meson04'bes'fah 



8 Frederick A. Harris 

as a square centered at (1.540, 1.540) GeV/c 2 , and we set an upper limit at the 90% 
confidence level (C.L.) on the branching ratio: 

B(ip(2S) -> 00 -> K%pK~n + K%pK + n) < 0.84 x 10~ 5 . 

Another possibility is that the ip(2S) decays to only one 8 or 9 state. To deter- 
mine the number of 0(1540) events from single or production, we count the 
number of events within regions of 1.52 - 1.56 GeV/c 2 in the projections of Fig.0 
and set upper limits, shown in Table |3] 



Tabic 3. Summary of upper limits. 



Decay mode 


i>(2S) 


J/i/> 


-> 66 -+ K^pK-fix + K^pK+n 


0.88 x 10" b 


1.1 x 10" 


b 


-> BK-n -> KlpK-fi 


1.0 x 10" 5 


2.1 x 10" 


-5 


-> OK+n -> K%pK+n 


2.6 x 10~ 5 


5.6 x 10" 


-5 


-» K% P e -> K° s pK-n 


0.60 x 10~ 5 


1.1 x 10" 


-5 


XgpO -> K° s pK+n 


0.70 x 10" 5 


1.6 x 10" 


-5 



For the decays of J/ip — > KgpK~n and KgpK + n, we use the same criteria 
and analysis method as those used for the ip{2S) data to study possible 0(1540) 
production. There is no significant 0(1540) signal, and we determine upper limits 
on the branching fractions at the 90% C.L., shown in Table El Full details may be 
found in Ref. El 

6. Enhancement in the pA invariant mass spectrum in 
J /ip — ► pK~A and in ip(2S) — ► pK~A decays 

An anomalous enhancement near the mass threshold in the pp invariant mass spec- 
trum was observed by the BES II experiment in J/ip — > jpp decays. ED This en- 
hancement can be fitted with an S-wave Breit-Wigner resonance function with a 
mass around 1860 MeV and a width T < 30 MeV, and has been interpreted as a 
possible baryonium state. ^^Similar pp mass-threshold enhancements have be en ob- 
served in the decays B + — > K + pp and B° — > D°pp by the Belle Collaboration. 
It is, therefore, of special interest to search for possible resonant structures in other 
baryon-antibaryon final states. The Belle Collaboration recently observed a near- 
threshold enhancement in the pA mass spectrum from B — > pAir decays. We 
report the observation of an enhancement near threshold in the pA invariant mass 
spectrum in J /if) — > pK~A and in ip(2S) — > pK~A decays. The results are based 
on an analysis of 5.8 x 10 7 J/ip and 1.4 x 10 7 tp(2S) decays detected in BESII. 

The J/ip — > pK~A candidate events are required to have four good charged 
tracks with total charge zero. Events are subjected to a four-constraint (4C) kine- 
matic fit with the corresponding mass assignments for each track. For events with 
ambiguous particle identification, all possible 4C combinations are formed, and the 
combination with the smallest x 2 is chosen. A sample of 5421 J/ip — ► pK~A candi- 
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dates survive the final selection. Monte Carlo studies indicate that the background 
in the selected event sample is at the 1 ~ 2% level. 



^200 




M pa " m P " m A (GeV) 



0.1 



COS0„ 



Fig. 7. (a) The points with error bars indicate the measured pA mass spectrum; the shaded 
histogram indicates phase space MC events (arbitrary normalization), (b) The Dalitz plot for the 
selected event sample, (c) A fit (solid line) to the data. The dotted curve indicates the Breit-Wigner 
signal and the dashed curve the phase space 'background.' (d) The cos6 p distribution under the 
enhancement, the points are data and the histogram is the MC (normalized to data) 



The pA invariant mass spectrum for the selected events is shown in Fig. Ufa), 
where an enhancement is evident near the mass threshold. No corresponding struc- 
ture is seen in a sample of J/ip — > pK~ A MC events generated with a uniform 
phase space distribution. The pK~A Dalitz plot is shown in Fig.EJb). In addition 
to bands for the well established A* (1520) and A* (1690), there is a significant N* 
band near the K~A mass threshold, and a pA mass enhancement, isolated from the 
A* and N* bands, in the right-upper part of the Dalitz plot. 

This enhancement can be fit with an acceptance weighted S-wave Brcit-Wigner 
function, together with a function describing the phase space contribution, as shown 
in Fig. |7Jc). The fit gives a peak mass of m — 2075 ± 12 ± 5 MeV and a width 
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r = 90±35±9 MeV. The significance is about 7<r. A P-wave Breit-Wigner resonance 
functions can also fit the enhancement. The cos# p distribution, shown in Fig.Efd), 
where 9 P is the decay angle of p in the pA CM frame, agrees well with that of a MC 
sample of J/ip — > KX — ► i^TpA. Since the MC cos8 p distribution is generated as a 
uniform S-wave distribution and the detected MC distribution agrees with data in 
Fig-EJd), the observed distribution for the enhancement is consistent with S-wave 
decays to pA. 

Evidence of a similar enhancement is observed in yj(2S) — > pK~ A when the 
same analysis is performed on the ip(2S) data sample. More detail can be found in 

Ref.na 
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